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I.  INTRODUCTION 


In  1965  a paper  appeared  in  the  Tenth  Symposium  (International 
on  Combustion)  by  J.  George  Sotter  of  Sheffield  University,  Sheffield, 
England  in  which  the  author  constructed  profiles  of  velocity,  and 
temperatures  vs^  distance  for  the  induction  and  explosion  zones  above 
cordite  burning  at  1000  psi  (7  MP4) . For  this  purpose  it  was  necessary 
to  select  pertinent  reactions  and  to  assemble  reaction  rate  expressions 
for  them.  The  author  selected  a total  of  twenty  reactions  and  carried 
out  computer  calculations  using  them. 

In  the  present  report  the  rate  constant  data  used  are  examined  in 
the  light  of  current  knowledge  of  these  reactions  and,  where  necessary, 
revised  values  are  suggested. 

II.  SOURCES  OF  REACTION  RATE  CONSTANTS 

It  is  an  unusual,  and  unfortunate  aspect  of  Sotter' s paper  that 
not  only  were  the  pertinent  rate  constant  data  not  tabulated,  but  the 
reference  source*-  was  a paper  given  at  a meeting  in  1962  of  the  Western 
States  Section  of  the  Combustion  Institute.  This  paper  was  not  pub- 
lished later  in  a journal  and  considerable  library  time  was  spent 
trying  to  find  it.  A copy  was  finally  obtained  through  the  kindness 
of  Ms.  Rebecca  Palmer,  Secretary  to  Professor  Richard  Sawyer,  Depart- 
ment of  Mechanical  Engineering,  University  of  California  who  had  pre- 
sided at  the  particular  meeting.  The  reactions  considered  by  Sotter 
and  the  rate  constants  he  used  are  shown  in  Table  I. 

Initially  it  was  planned  to  survey  the  literature  for  the  indi- 
vidual reactions  used  by  Sotter  and  to  evaluate  the  results  of  this 
survey.  It  soon  became  evident  that  this  was  a task  of  very  consider- 
able magnitude.  Fortunately  the  examination  of  the  literature  re- 
vealed the  existence  of  critical  surveys  of  reaction  rate  data  which 
included  most  of  the  reactions  of  interest.  Since  the  surveys  origi- 
nated in  laboratories  where  considerable  expertise  existed  in  connection 
with  reactions  of  the  type  of  interest  here,  the  values  reported  here 
are  for  the  most  part  taken  from  those  surveys.  In  Table  I the  values 
recommended  are  shown  along  with  those  used  by  Sotter. 


J.  G.  Sotter,  Tenth  Symposium  (International)  on  Combustion,  1405, 
(1965). 

G.  E.  Ellis  and  G.  S.  Bahn,  "Engineering  Selection  of  Reaction 
Rate  Constants  for  Gaseous  Chemical  Species  at  High  Temperatures ", 
1962  Fall  Meeting,  Western  States  Section,  The  Combustion  Institute. 
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III.  REMARKS  ON  CURRENT  CONSTANTS  AND  THOSE  USED  BY  SOTTER 


Of  the  twenty  reactions  cited  by  Sotter,  only  seventeen  are  inde- 
pendent. Reactions  (4b),  (5b),  and  (6b)  of  Table  I were  assigned  con- 
stants on  the  basis  of  applying  an  arbitrary  "catalyst-efficiency" 
factor  to  the  rate  constants  for  reactions  (4) , (5) , and  (6) , 
respectively. 

In  considering  such  a large  number  of  reactions,  it  is  a problem 
deciding  where  to  start.  It  is  an  additional  problem  to  decide  how 
significant  discrepancies  between  current  values  and  those  used  by 
Sotter  are.  In  his  paper  Sotter  emphasized  the  importance  of  the  "water- 
gas  equilibrium"  reactions,  reactions  (1)  and  (7)  of  Table  I and  of 
reaction  (17) , the  termolecular  reaction  between  nitric  oxide  and  hydro- 
gen. He  states,  in  addition,  that  reactions  (8),  (9),  (10),  (12),  and 
(14)  could  probably  be  omitted  without  loss  of  accuracy. 

It  has  therefore  seemed  desirable  to  add  a few  comments  about 
these  reactions.  Since  no  new  data  have  been  found  for  reaction  (17), 
comments  will  be  made  only  for  the  other  two. 

Reaction  (1) 

H + H^O  ■+  H2  + OH 

3 4 

has  been  considered  by  the  Leeds  group  and  by  Wilson  . The  value 
recommended  by  the  former  for  the  reaction  for  which  Sotter  lists  the 
value,  i.e. , the  back  reaction 

OH  + H2  -*■  H + H.,0 

13  3 -1  -1 

is  k = 2.19  x 10  exp  (-5150/RT)  cm  mole  sec  *35%  and  is  given  for 
the  temperature  range  300-3000K.  Wilson4  recommends 

13  3 -1  -1 

k » 2.3  x 10  exp(-5200/RT)  cm  mole  sec  for  the  temperature 

range  300-2000K.  The  two  values  are  thus  in  close  agreement  and  differ 
somewhat  from  the  value  used  by  Sotter.  At  1000K  for  example  the  ex- 
pression used  by  Sotter  yields  a rate  constant  some  thirty  times  as 
large  as  that  derived  from  the  above  values. 


(a)  D.  L.  Baulch  et  al. , "High  Temperature  Reaction  Rate  Data ", 

No.  1 , Dept  of  Physical  Chemistry , The  University , Leeds , 2,  England ; 
May  1968  (b)  ibid , No.  2;  Nov  1968  (c)  ibid  No.  3 ; Apr  1969  (d)  ibid , 
No.  4;  Dec  1969  (e)  ibid , No.  S;  Jul  1970. 

^ W.  E.  Wilson , Jr. , "A  Critical  Review  of  the  Gas  Phase  Reaction 

Kinetics  of  the  Hydroxyl  Radical ",  J.  Phys.  Chem.  Ref.  Data , 11 , S35 

(1972). 
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Perhaps  the  most  serious  revision  of  all  occurs  for  reaction  (7) 

CO  + OH  -*•  C02  + H 

The  value  used  by  Sotter  was  taken  from  the  work  of  Westenberg  and 
Favin^  at  the  Applied  Physics  Laboratory, 

k^  = 10*^exp(-10,000/RT)  cm^mole~*sec”* 

This  expression  was  based  on  measurements  made  in  the  authors'  labora- 
tory at  elevated  temperatures  and  on  data  reported  for  room  temperature 
measurements  in  the  Russian  literature^.  Further  studies  in  Westen- 
berg’s  laboratory?  showed  the  Russian  work  to  be  in  serious  error. 
Reaction  rate  constants  have  been  evaluated  as  shown  below. 

Measured  Rate  Constant  for  CO  + OH 

3.-1  -1 

in  cm  mole  sec 


T(K)  ^f 


298 

00 

o 

X 

10*° 

523 

8.7 

X 

io10 

915 

13.1 

X 

io10 

The  data  indicate  an  activation  energy  of  zero.  A general  expression 
of  the  following  type  can  be  written 

13 

. _ 2.5  x 10 

1 " Tn[l-exp(-hvi/kT)] 

where  five  vibrations  must  be  considered  for  v. : 1670  cm  * (C-0  stretch), 

1300  cm-l  (0-H  bend),  830  cm-1  (H  out-of-plane1bend) . Values  are  given 
in  the  reference  for  the  denominator  from  300-1000K.  Wilson4  recommends 
for  the  range  300-2000K 

k = 3.1  x 10**  exp-(300/T)cm3mole~*sec~* 


A.  A.  Westenberg  and  S.  Favin,  Ninth  Symposium  (International)  on 
Combustion,  Academia  Press , 1963,  p.  78S. 

L.  I.  Avramenko  and  R.  V.  Lorentso , Zhur.  Fiz.  Khim . , 2£,  207  (1950). 
A.  A.  Westenberg  and  N.  DeHaas,  J.  Chem.  Phys 58j  406  (1973). 
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which  yields  k = 2.3  x loj^sec”^  at  1000K.  Westenberg  and  DeHaas 
report  k = 1.4  x lO^sec"  . The  Leeds  Group  3a  recommends 

k = 5.6  + .8  x 10nexp(-1080  * SOO/RTjcnwnole^sec"1 

At  1000K  this  yields  k = 3.4  x lO^sec-1  which  is  also  in  good  agree- 
ment. There  appears  to  be  little  to  choose  from  among  these.  I favor 
the  data  of  Westenberg  and  DeHaas  but  would  like  to  point  out  that  the 
other  expressions  yield  rather  similar  results. 

Where  the  remaining  rate  constants  show  differences,  they  are 
much  smaller  and  I am  not  able  to  evaluate  their  importance  for  the 
results  obtained  by  Sotter.  The  groups  cited  seem  most  reliable,  and 
since  no  serious  differences  exist  for  those  rate  constants  for  which 
values  are  cited  independently  by  both  of  them,  the  recommended  values 
should  be  reliable. 

Note  that  in  Table  I rate  constants  are  cited  for  either  the  for- 
ward or  back  reaction  in  the  reversible  reactions  but  not  for  both. 
Sotter  obtained  equilibrium  constants^  for  all  reactions  except  No.  7 
of  Table  I.  The  equilibrium  constants  in  reference  2 were  calgulated 
from  data  tabulated  in  the  JANAF  Interim  Thermochemical  Tables  . 
Equilibrium  constants  are  generally  available  in  more  reliable  form 
than  are  rate  constants  and,  since  the  cited  source  is  one  whose 
authority  is  widely  accepted,  I have  not  felt  it  necessary  or  desirable 
to  calculate  and  tabulate  either  equilibrium  constants  or  the  other 
rate  constant  for  each  reversible  reaction.  It  is  clear,  of  course, 
that  the  combination  of  an  incorrect  k^.  and  a correct  equilibrium  con- 
stant leads  to  an  incorrect  k,  . 


D.  R.  Stull  et  als  JANAF  Interim  Thermochemical  Tables , Thermal 
Laboratory , Dow  Chemical  Co.,  Midland , MI,  Dec  31,  1960,  eq  seq. 
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Table  I.  Reaction  Rates  Constants  Used  by  Sotter  and 
Values  Recommended  from  More  Recent  Compilations 

(Values  used  by  Sotter  were  taken  from  reference  2 except  for  Reaction  (7)  which  was  taken  from 
reference  5,  and  Reaction  (17)  which  is  given  by  Sotter  along  with  the  reference  source.  For 
the  Recommended  Rate  Constants,  all  values  were  taken  from  reference  3 except  for  Reaction  (7), 
which  is  discussed  earlier  in  this  report,  and  Reaction  (13)  for  which  the  source  is  reference  9. 

No.  Reaction  Sotter's  Value  Current  Recommendation 


(1) 

H + 

H20 

H9  + OH 

(2) 

H + 

°2 

0 + OH 

(3) 

0 + 

H2 

■+  OH  + H 

(4) 

h2o 

+ M 

-*•  H + OH  + M 

(4b) 

H + 

OH 

+ NO  H^O  + NO 

(5) 

H2  ■ 

* M 

-*■  211  + M 

(5b) 

H + 

H + 

NO  -*■  H2  + NO 

(6) 

°2  + 

M - 

20  + M 

(6b) 

0 + 0 + 

NO  - 0o  + NO 

(7) 

CO  + 

OH 

->  C02  + II 

(8) 

NO  + 

M -*• 

N + 0 + M 

(9) 

N2  ♦ 

0 -► 

NO  + N 

kb=3xl014exp(-2500/RT) 

kf=4xl014exp(-18000/RT) 

kf=3xl014exp(-6000/RT) 

kf-1023T'1*5exp(-114730/RT) 

k^  for  (4)  multiplied  by  factor 

kf=1021T"1,5exp(-103240/RT) 

k^  for  (5)  multiplied  by  factor 

kf=1.1369xl025T'2,5exp(-118000/RT) 

k^  for  (6)  multiplied  by  factor 

kf=1013exp(-10000/RT) 

kf=5. 18x1021T_1 •5exp(-150000/RT) 
kf=5.5xl0nT°'5exp(-72000/RT) 


kb=2.2xl013exp(-5150  nSO/RT) 
kt=2.2xl014± -3exp(-16800  ±500/RT) 
kf=1.7xl013±,2exp(-9450  *250/RT) 
akf=l .6xl018exp(-119000/RT) 
no  new  data 
no  new  data 
no  new  data 
no  new  data 
no  new  data 

13 

2.5x10 

k Tit  [l-exp(-hv  i/kT)) 
kf=5.6*.8xl0nexp(-1080*500/RT) 
kf=l.lxl021T-1,5exp(-150000/RT) 
Ckf=1.36xl014exp(-75000/RT) 


No. 

Reaction 

Sotter's  Value 

Current  Recommendation 

(10) 

NO  + 0 -*•  0?  + N 

kf=3.2xl09exp(-39100/RT) 

kf=1.5xl09Texp(-38640/RT) 

(11) 

N?  + 0o  ->  2N0 

kf=9.1xl024T'2,5exp(-128500/RT)no  new  data 

(12) 

0 + N,0  2N0 

kb=1.0xlonexp(-15500/RT) 

kb=2.3xl013exp(-24100/RT) 

(13) 

N20  + M N2  + 0 + M 

kf=2.07xl022T‘2exp(-62000/RT) 

kf=4.88xl014exp(-58000/RT) 

(14) 

N?  + °2  - M20  + 0 

kf =4 . 0xl09exp ( -93200/RT) 

no  new  data 

(15) 

0 + H„0  20H 

kb=3xl014exp(-2S00/RT) 

kb=5.8xl012*  *45exp(-780±500/RT) 

(16) 

OH  + M ->  H + 0 + M 

kf=1021T'1,5exp(-101000/RT) 

no  new  data 

(17) 

2N0  + H2  - 2HN0  (+H2 

1 ft 

kf=5xl0  exp(-47000/RT) 

no  new  data 

- 2H20  + N2) 

a.  5 3 

Estimated  from  an  estimate  in  reference  3 that  E = 119,000  and  that  at  200K,  k=3.4xl0  cm  mole 

^v=cf.  remarks  in  text 

c 

Accuracy  given  as±100% 


K.  Schofield,  "Evaluated  Chemical  Kinetic  Rate  Constants  for  Various  Gas  Phase  Reactions", 
J.  Phys.  Chem.  Ref.  Data.,  2,  25  (1973). 
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